- Author:
- Weiwei Ai <wai484@aucklanduni.ac.nz>
- Date:
- 2024-04-18 18:08:11+12:00
- Desc:
- Add BG for steady-state simulation (i.e, a longer pulse)
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/b65/rawfile/c8bc4b889b3850d67cea87048416591cdb681ffd/Electrogenic cotransporter/SGLT1_BG.txt
def model SGLT1_BG as
def import using "./units.cellml" for
unit J_per_K_mol using unit J_per_K_mol;
unit J_per_mol using unit J_per_mol;
unit fA using unit fA;
unit fmol using unit fmol;
unit per_fmol using unit per_fmol;
unit C_per_mol using unit C_per_mol;
unit fmol_per_s using unit fmol_per_s;
enddef;
def import using "params_BG.cellml" for
comp params_BG using comp params_BG;
enddef;
def comp SGLT1_BG as
var t: second;
var R: J_per_K_mol{ pub: in};
var T: kelvin{ pub: in};
var F: C_per_mol{ pub: in};
var K_Nai: per_fmol{ pub: in};
var q_Nai: fmol{ pub: in};
var K_Nao: per_fmol{ pub: in};
var q_Nao: fmol{ pub: in};
var K_Glci: per_fmol{ pub: in};
var q_Glci: fmol{ pub: in};
var K_Glco: per_fmol{ pub: in};
var q_Glco: fmol{ pub: in};
var K_1: per_fmol{ pub: in};
var q_init_1: fmol{ pub: in};
var K_2: per_fmol{ pub: in};
var q_init_2: fmol{ pub: in};
var K_3: per_fmol{ pub: in};
var q_init_3: fmol{ pub: in};
var K_4: per_fmol{ pub: in};
var q_init_4: fmol{ pub: in};
var K_5: per_fmol{ pub: in};
var q_init_5: fmol{ pub: in};
var K_6: per_fmol{ pub: in};
var q_init_6: fmol{ pub: in};
var V0_Vm: volt{ pub: in};
var z_zf1: dimensionless{ pub: in};
var z_zr1: dimensionless{ pub: in};
var z_zf6: dimensionless{ pub: in};
var z_zr6: dimensionless{ pub: in};
var kappa_r1: fmol_per_s{ pub: in};
var kappa_r2: fmol_per_s{ pub: in};
var kappa_r3: fmol_per_s{ pub: in};
var kappa_r4: fmol_per_s{ pub: in};
var kappa_r5: fmol_per_s{ pub: in};
var kappa_r6: fmol_per_s{ pub: in};
var kappa_r7: fmol_per_s{ pub: in};
var mu_Nai: J_per_mol;
var mu_Nao: J_per_mol;
var mu_Glci: J_per_mol;
var mu_Glco: J_per_mol;
var mu_1: J_per_mol;
var v_1: fmol_per_s;
var mu_2: J_per_mol;
var v_2: fmol_per_s;
var mu_3: J_per_mol;
var v_3: fmol_per_s;
var mu_4: J_per_mol;
var v_4: fmol_per_s;
var mu_5: J_per_mol;
var v_5: fmol_per_s;
var mu_6: J_per_mol;
var v_6: fmol_per_s;
var V_Vm: volt;
var V_zf1: volt;
var I_zf1: fA;
var mu_zf1: J_per_mol;
var v_zf1: fmol_per_s;
var V_zr1: volt;
var I_zr1: fA;
var mu_zr1: J_per_mol;
var v_zr1: fmol_per_s;
var V_zf6: volt;
var I_zf6: fA;
var mu_zf6: J_per_mol;
var v_zf6: fmol_per_s;
var V_zr6: volt;
var I_zr6: fA;
var mu_zr6: J_per_mol;
var v_zr6: fmol_per_s;
var A_f_r1: J_per_mol;
var v_r1: fmol_per_s;
var A_r_r1: J_per_mol;
var A_f_r2: J_per_mol;
var v_r2: fmol_per_s;
var A_r_r2: J_per_mol;
var A_f_r3: J_per_mol;
var v_r3: fmol_per_s;
var A_r_r3: J_per_mol;
var A_f_r4: J_per_mol;
var v_r4: fmol_per_s;
var A_r_r4: J_per_mol;
var A_f_r5: J_per_mol;
var v_r5: fmol_per_s;
var A_r_r5: J_per_mol;
var A_f_r6: J_per_mol;
var v_r6: fmol_per_s;
var A_r_r6: J_per_mol;
var A_f_r7: J_per_mol;
var v_r7: fmol_per_s;
var A_r_r7: J_per_mol;
var q_1: fmol{ init: q_init_1};
var q_2: fmol{ init: q_init_2};
var q_3: fmol{ init: q_init_3};
var q_4: fmol{ init: q_init_4};
var q_5: fmol{ init: q_init_5};
var q_6: fmol{ init: q_init_6};
mu_Nai = R*T*ln(K_Nai*q_Nai);
mu_Nao = R*T*ln(K_Nao*q_Nao);
mu_Glci = R*T*ln(K_Glci*q_Glci);
mu_Glco = R*T*ln(K_Glco*q_Glco);
mu_1 = R*T*ln(K_1*q_1);
ode(q_1,t) = v_1;
mu_2 = R*T*ln(K_2*q_2);
ode(q_2,t) = v_2;
mu_3 = R*T*ln(K_3*q_3);
ode(q_3,t) = v_3;
mu_4 = R*T*ln(K_4*q_4);
ode(q_4,t) = v_4;
mu_5 = R*T*ln(K_5*q_5);
ode(q_5,t) = v_5;
mu_6 = R*T*ln(K_6*q_6);
ode(q_6,t) = v_6;
V_Vm=V0_Vm;
mu_zf1=z_zf1*F*V_zf1;
I_zf1=z_zf1*F*v_zf1;
mu_zr1=z_zr1*F*V_zr1;
I_zr1=z_zr1*F*v_zr1;
mu_zf6=z_zf6*F*V_zf6;
I_zf6=z_zf6*F*v_zf6;
mu_zr6=z_zr6*F*V_zr6;
I_zr6=z_zr6*F*v_zr6;
v_r1 = kappa_r1*(exp(A_f_r1/(R*T)) - exp(A_r_r1/(R*T)));
v_r2 = kappa_r2*(exp(A_f_r2/(R*T)) - exp(A_r_r2/(R*T)));
v_r3 = kappa_r3*(exp(A_f_r3/(R*T)) - exp(A_r_r3/(R*T)));
v_r4 = kappa_r4*(exp(A_f_r4/(R*T)) - exp(A_r_r4/(R*T)));
v_r5 = kappa_r5*(exp(A_f_r5/(R*T)) - exp(A_r_r5/(R*T)));
v_r6 = kappa_r6*(exp(A_f_r6/(R*T)) - exp(A_r_r6/(R*T)));
v_r7 = kappa_r7*(exp(A_f_r7/(R*T)) - exp(A_r_r7/(R*T)));
v_1 = -v_r1 +v_r6;
v_2 = v_r1 -v_r2 -v_r7;
v_3 = v_r2 -v_r3;
v_4 = v_r3 -v_r4;
v_5 = v_r4 -v_r5 +v_r7;
v_6 = v_r5 -v_r6;
v_zf1 = v_r1;
V_zf1 = V_Vm;
v_zr1 = v_r1;
V_zr1 = V_Vm;
v_zf6 = v_r6;
V_zf6 = V_Vm;
v_zr6 = v_r6;
V_zr6 = V_Vm;
A_f_r1 = 2{dimensionless}*mu_Nao +mu_1 -mu_zf1;
A_r_r1 = mu_2 +mu_zr1;
A_f_r2 = mu_Glco +mu_2;
A_r_r2 = mu_3;
A_f_r3 = mu_3;
A_r_r3 = mu_4;
A_f_r4 = mu_4;
A_r_r4 = mu_Glci +mu_5;
A_f_r5 = mu_5;
A_r_r5 = 2{dimensionless}*mu_Nai +mu_6;
A_f_r6 = mu_6 -mu_zf6;
A_r_r6 = mu_1 +mu_zr6;
A_f_r7 = mu_2;
A_r_r7 = mu_5;
enddef;
def map between params_BG and SGLT1_BG for
vars R and R;
vars T and T;
vars F and F;
vars K_Nai and K_Nai;
vars q_Nai and q_Nai;
vars K_Nao and K_Nao;
vars q_Nao and q_Nao;
vars K_Glci and K_Glci;
vars q_Glci and q_Glci;
vars K_Glco and K_Glco;
vars q_Glco and q_Glco;
vars K_1 and K_1;
vars q_init_1 and q_init_1;
vars K_2 and K_2;
vars q_init_2 and q_init_2;
vars K_3 and K_3;
vars q_init_3 and q_init_3;
vars K_4 and K_4;
vars q_init_4 and q_init_4;
vars K_5 and K_5;
vars q_init_5 and q_init_5;
vars K_6 and K_6;
vars q_init_6 and q_init_6;
vars V0_Vm and V0_Vm;
vars z_zf1 and z_zf1;
vars z_zr1 and z_zr1;
vars z_zf6 and z_zf6;
vars z_zr6 and z_zr6;
vars kappa_r1 and kappa_r1;
vars kappa_r2 and kappa_r2;
vars kappa_r3 and kappa_r3;
vars kappa_r4 and kappa_r4;
vars kappa_r5 and kappa_r5;
vars kappa_r6 and kappa_r6;
vars kappa_r7 and kappa_r7;
enddef;
enddef;
